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A fundamental problem in plasma and astrophysics is the interaction between energetic particles 
and magnetized plasmas. In the current paper we focus on particle diffusion across the guide 
magnetic field. It is shown that the perpendicular diffusion coefficient depends only on the parallel 
diffusion coefficient and the Kubo number. Therefore, one can find four asymptotic limits depending 
on the values of these two parameters. These regimes are the quasilinear limit, the Kadomtsev & 
Pogutse limit, the scaling of Rechester & Rosenbluth, and the scaling found by Zybin & Istomin. In 
the current article we focus on the Rechester & Rosenbluth scenario because this was not discovered 
before in the context of collisionless plasmas. Examples and applications are discussed as well. We 
show that an energy independent ratio of perpendicular and parallel diffusion coefficients can be 
found and that this ratio can be very small but also close to unity. This is exactly what one observes 
in the solar wind. 

PACS numbers: 47.27.tb, 96.50.Ci, 96.50.Bh 


I. INTRODUCTION 

A fundamental problem in physics is the stochastic 
motion of particles. Originally the stochastic motion of 
Brownian particles was investigated by R. Brown and A. 
Einstein (see Refs. [1] and [2]). Later S. Chandrasekhar 
presented a more general discussion of stochastic prob¬ 
lems in Physics and Astronomy (see Ref. [3]). In the re¬ 
cent decades the interaction between energetic particles 
and a turbulent magnetized plasmas became a topic of 
great interest (see, e.g., Ref. [4]). An interesting example 
is the motion of alpha particles or accelerated ion beams 
in a Tokamak (see, e.g., Ref. [5]). In space science and 
astrophysics, cosmic rays propagate through the plasma 
of the interstellar space or the solar wind. In their report 
to the National Academy of Science , Turner et al. (see 
Ref. [6]) listed the eleven fundamental science questions 
for the 21st century. One of these questions is the nature 
of cosmic ray sources. One ingredient to the solution 
of this puzzle is the understanding of energetic particle 
diffusion because this plays a crucial role in shock accel¬ 
eration scenarios. 

If charged particles move through the plasma, they ex¬ 
perience scattering due to the interaction with turbulent 
magnetic fields. In addition to such stochastic magnetic 
fields, in the following represented by the symbol SB, one 
also finds a non-vanishing mean (guide) magnetic field 
Bq = Bae z . Since the latter field breaks the symmetry 
of the physical system, one has to distinguish between 
diffusion of particles along and across the mean field. In 
analytical treatments of the transport, it appears that 
perpendicular diffusion is more difficult to understand 
(see, e.g., Ref. [7]). 

Recently the so-called Unified Non-Linear Transport 
(UNLT) theory was presented in Ref. [8]. The valid¬ 
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ity of the latter theory was tested by comparing its re¬ 
sults with test-particle simulations. Remarkable agree¬ 
ment was found for a variety of turbulence models such 
as the slab/2D model and nosiy reduced magnetohydro¬ 
dynamic (MHD) turbulence (see Refs. [9] and [10]). In 
Ref. [11] it was shown that UNLT theory agrees also with 
test-particle simulations performed for models based on 
Goldreich & Sridhar scaling (see Ref. [12]). 

In the current paper we show that the perpendicu¬ 
lar diffusion coefficient depends only on two parameters, 
namely the parallel mean free path and the Kubo num¬ 
ber. Furthermore, a new formula for the perpendicular 
diffusion coefficient is derived from UNLT theory corre¬ 
sponding to the scaling obtained by Rechester & Rosen¬ 
bluth (see Refs. [13] and [14]). After discussing examples, 
it is argued that the latter transport regime is important 
for explaining solar wind observations. 


II. PERPENDICULAR DIFFUSION 

The UNLT theory provides the following nonlinear in¬ 
tegral equation for the perpendicular diffusion coefficient 

_ a 2 v 2 f _ P X x(k) _ 

3 B% J F(k\\,k±) + (4/3)Kj_fcj_ + t>/A|| 

with F(k\\,kj_) = (r' 2 fcj|)/(3ftj_fc^). Here we have 
used the parallel and perpendicular components of the 
wave vector k\\ and k±, the magnetic correlation tensor 

Pmn{k) = (• 5B m (k)SB*(k )), the parallel diffusion coef¬ 
ficient of the particle K||, the parallel mean free path 
Ay = 3k||/u, and the particle speed v. We have also 
used the order one parameter a 2 as in Ref. [15]. Eq. (1) 
is valid for axis-symmetric static turbulence, SB Z -C Bq, 
and a constant guide field. 

In the following we assume that the combined correla- 
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tion function P(k) = P xx (k) + P yy (k ) has the form 

p{k)=l ll llSB 2 x f(k ll l ll ,k ± l ± ) (2) 

with the dimensionless function f(x,y ) which depends 
only on x = fc||Z|| and y = k±l±. Furthermore, this 
function decays with increasing parallel and perpendicu¬ 
lar wavenumbers. Therefore, 1 1 | and l± are characteris¬ 
tic length scales for the decorrelation of the turbulence. 
The form (2) was chosen so that f d 3 k P{k) = SB 2 = 
SB 2 + SBy. Therefore, the function / satisfies the nor¬ 
malization condition 1 = 27r dy J^° dx yf(x,y). To 
proceed we use the diffusion ratio D := = 

(Aj_Zj?)/(A||Zj_) and the Kubo number 


K = h 

l± Bq 


( 3 ) 


Therewith, Eq. (1) can be written as 

nOO />00 

D = 2tto, 2 K 2 / dy dxyf(x,y) 
Jo Jo 


x 2 /y 2 + D + (2DA||) 2 / (3Z,,) 2 y 2 ’ 

Obviously, there are only two parameters controlling the 
diffusion ratio D , namely the parallel mean free path nor¬ 
malized to the parallel scale An/Z|| and the Kubo number 
K. We conclude that there are four different asymptotic 
limits and therewith four different transport regimes. 
Those are discussed in the following. All limits are com¬ 
pared with each other in Table I. 


A. The quasilinear regime 


First we consider the limit An/Zm —> oo corresponding 
to the case that pitch-angle scattering and therewith par¬ 
allel diffusion are suppressed. In this case it is convenient 
to go back to Eq. (1) and set v/X\\ = 0. In this limit Eq. 
(1) has the solution = vkfl /2 where the field line 
diffusion coefficient kfl is given by 


k-fl = 


B 2 

n o 


d 3 k P xx (kj 


KFLkj_ 


kn + (kfzA’1) 


( 5 ) 


in agreement with the theory developed by Matthaeus et 
al. (see Ref. [16]). Here we have introduced the effective 
guide field B 0 = B 0 /a. By using again the form (2) 
and by introducing the dimensionless field line diffusion 
coefficient D = KflI\\/1±, Eq. (5) becomes 


D = 2t tK 2 




dx yf(x,y) 



( 6 ) 


We find that D depends only on one single parameter 
and that is the Kubo number K. One can obtain two 


different asymptotic solutions of the latter integral equa¬ 
tion by considering small and large Kubo numbers. In 
the limit K ->0we expect D —> 0. By using the relation 
lim^o £,/(x 2 + £ 2 ) = nS(x), with the Dirac delta S(x), 
Eq. (6) becomes 


/»oo 

D = 2tt 2 K 2 / dy yf(x = 0, y) (7) 

Jo 

and the parameter D , and therewith the field line diffu¬ 
sion coefficient kfl, scales with K 2 corresponding to the 
well-known quasilinear scaling (see, e.g., Ref. [17]). 


B. The Kadomtsev & Pogutse scaling 

In the limit of large Kubo numbers K —> oo we expect 
D —> oo. In this case Eq. (6) becomes 

/»00 f‘ OO 

D 2 =2ttK 2 dy / dx y~ 1 f(x,y) (8) 

Jo Jo 

and we find the scaling D oc K. The latter limit is 
sometimes called the nonlinear regime or the Bohrn limit 
of field line diffusion and was originally obtained by 
Kadomtsev & Pogutse (see Ref. [18]). 


C. The Zybin &: Istomin scaling 

In the limit A||/Z|| —> 0, Eq. (4) becomes 

r 00 r°° n 

D = 2ira 2 K 2 / dy / dx yf{x,y) 9 . (9) 

Jo Jo x 2 /y 2 + D 

For large Kubo numbers we expect that the diffusion pa¬ 
rameter D becomes large as well and Eq. (9) simplifies 
to 

/>00 nOO 

D = 2na 2 K 2 / dy dx yf(x,y). (10) 
Jo Jo 

Obviously D oc K 2 . Eq. (10) can be rewritten as k± = 
a 2 K\\SB 2 /Bq- The latter result is also known in diffusion 
theory (see Refs. [19] and [7]) and we refer to it as the 
Zybin & Istomin scaling. 


D. The collisionless Rechester &; Rosenbluth 
scaling 

For small Kubo numbers we use the relations 
limc^o i/{x 2 + £ 2 ) = irS(x) and 5 (az) = 8{z)/\a\ to 
derive from Eq. (9) 

VD = 2n 2 a 2 K 2 / dy y 2 f(x = 0,y). (11) 

Jo 

Assuming that the y -integral is convergent, we find the 
scaling D oc A' 4 . To proceed we define the dimensionless 
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parameter y 2 := 2i t 2 dy y 2 f(x = 0 ,y). Furthermore, 
we replace the parameters D and K by standard quanti¬ 
ties to obtain the scaling 


k± 

K \\ 


2 4 
= 7 a 


I 5Bj 


B 4 0 


( 12 ) 


The latter limit was derived from the UNLT theory. In 
order to understand this result, we compare our find¬ 
ing with the formulas discussed in Ref. [13]. Accord¬ 
ing to the latter reference, Reclrester & Rosenbluth pro¬ 
posed k _l oc (/||K||dl3 2 )/(L/fRg) with the Kolmogorov 
length Lk■ If the latter parameter is replaced by Lk = 
(y/^l^B 2 ) / (4y/nl\\SB 2 ) (see, e.g., Ref. [20]), one obtains 


K± 


2\Z2tt 


l p_Bt 

B 4 


(13) 


which agrees with Eq. (12). By replacing D and K by 
standard parameters, Eq. (11) can be written as 



d 3 k P xx (j?J k±6(k\\). 


(14) 


The latter formula can easily be combined with any tur¬ 
bulence model as long as the occurring integrals are con¬ 
vergent. We refer to this limit as the collisionless Roch¬ 
ester & Rosenbluth (CLRR) scaling. 


III. ILLUSTRATIVE EXAMPLES 

A simple model for turbulence is the slab model where 
all wave vectors are assumed to be parallel wrt the mean 
field. It is well-known that perpendicular transport is 
subdiffusive in this case (see, e.g., Ref. [21]). Of course 
this subdiffusive behavior does only occur if the parallel 
mean free path is finite. In the asymptotic limit A|| —> oo, 
the usual quasilinear scaling can be obtained as discussed 
above. For a finite parallel mean free path, we find a 
transport process which is usually called compound dif¬ 
fusion (see Ref. [22] for details). In the following it 
is shown that diffusion is restored if the slab model is 
broadened. The broadening of a turbulence model with 
reduced dimensionality was discussed before in Ref. [23]. 
Physically this corresponds to the case where wave vec¬ 
tors are mainly oriented parallel wrt the mean field but 
weak fluctuations are taken into account. We refer to 
this model as the noisy slab model and we define it via 

p {k) = ^g(k l{ )e(l-k ± l ± ) (15) 

where we have used the Heaviside step function 0(;r). 
For the spectrum of the slab modes we employ g(k\\) = 

C{s)/{2k)SB\ [l + (fc|^||) 2 j~ s/2 (see, e.g., Ref. [24]). 
The normalization function C(s) = r(s/2)/[2y / 7rT((s — 
l)/2)] depends on the inertial range spectral index s and 


Gamma functions. By combining the noisy slab model 
with Eq. (14), one obtains 


K± 



l\\ SB 2 ' 

'll 


2 


(16) 


We can recover the case of pure slab turbulence by setting 
l± = oo and, thus, we find = 0 corresponding to 
subdiffusion. To estimate a number for the ratio k±/kh, 
we use 8B 2 /Bq « 0.5, l\\/l± ~ 0.75, and a 2 « 1. For 
s = 5/3 (see Ref. [25]) we find C(s = 5/3) « 0.12 and, 
thus, Eq. (16) provides k±/k.i\ « 0.005. 

As a second example we employ the Gaussian decorre¬ 
lation model used in Ref. [20] which is given by 


p(k) = 


(17) 


A problem of this model is the rapid decorrelation if the 
wave numbers are larger than the inverse scales. This 
behavior is not realistic in astrophysical scenarios such 
as the solar wind. In reality one expects a power law 
behavior with k~ 5 / 3 (see Ref. [25]). Such small scales, 
usually referred to as inertial range, and even smaller 
scales, are essential if it comes to pitch-angle scattering 
and therewith parallel diffusion (see, e.g., Refs. [24] and 
[7]). In the context of fusion plasmas, however, the form 
(17) is commonly used. Combining formula (14) with the 
spectrum (17), one finds 


= jj 5B i 
«|| V 4 ) l\ Bt ■ 


(18) 


The latter result is in agreement with the form (12) de¬ 
rived above. Furthermore, we can compare our find¬ 
ing with Eq. (13). We find agreement if we set a 2 = 
7r - 3 / 4 2 11 /4/3 po 0.95 confirming the assumption that a 2 
is indeed an order one constant. 


IV. SUMMARY AND CONCLUSION 

We have shown that the perpendicular diffusion co¬ 
efficient of energetic particles interacting with turbulent 
magnetic Helds depends only on two parameters, namely 
the ratio A||/Z|| and the Kubo number K. By considering 
asymptotic limits, we can identify four transport regimes, 
namely the quasilinear regime, the Kadomtsev & Pogutse 
limit, the Zybin & Istomin scaling, and the collisionless 
Rechester & Rosenbluth (CLRR) scaling. All cases are 
summarized in Table I. There we also list the corre¬ 
sponding formulas for the perpendicular diffusion coeffi¬ 
cient and the scaling behavior. 

In our work we focus on the CLRR scaling which was 
derived the first time from UNLT theory. This transport 
regime is obtained for short parallel mean free paths and 
small Kubo numbers. In this case the ratio k±/k\\ is con¬ 
stant and very small. Requiring convergence of the corre¬ 
sponding wavenumber integral, the perpendicular diffu¬ 
sion coefficient is in this case given by Eq. (14). We like 
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TABLE I. Transport regimes described by UNLT theory for extreme values of the parallel mean free path wrt the parallel scale 
A|| /Z|| and the Kubo number K. For each case we have listed the formula for the perpendicular diffusion coefficient k ± and the 
corresponding scaling behavior. All used parameters are explained in the main part of the text. 


Transport Regime 

Limit 

Perpendicular Diffusion Coefficient 

Scaling 

Quasilinear Diffusion 

Kadomtsev & Pogutse 

Rechester & Rosenblutli 

Zybin & Istomin 

A||//|| —y oo, Is. — y 0 

A||//|| —y oo, Is. — y oo 

A||//|| —y 0, Is. — y 0 

A||//|| —y 0, K — y oo 

«L=f^ f d 3 k P xx (k)5(k ||) 

°r - -i 1/2 

K± = j§- 0 [jd 3 kP xx (k)kf \ 

k± = [fd 3 k P xx (k)k ± 5(k l{ )] 2 

«l = f d 3 k P xx (k ) = a 2 tC||-g# 

\ 7 

A_l OC l M-5# 

II P>Q 

A ±( xA ,’Pfi 

II L ± 

\ \ SBl 

A_l oc 

II P>Q 

Compound Diffusion 

A||//|| finite, K = 0 

H (( Aa; )) 2 = KFLy/K\\/(Trt) 

Subdiffusive 


to point out that a collisionless adaption of Rechester & 
Rosenblutli was already presented in Ref. [26]. The letter 
result, however, disagrees with the formula derived in the 
current article and simulations of test-particle transport 
(see, e.g., [15]). 

A characteristic aspect of the quasilinear and the 
Kadomtsev & Pogutse regimes is the fact that the perpen¬ 
dicular diffusion coefficient does not depend on the mag¬ 
netic rigidity whereas in the other two cases, the ratio 
k_l/k|| is energy independent and the Kubo number con¬ 
trols the magnitude of the ratio k_l/k||. For large Kubo 
numbers k±/k\\ is close to unity as long as SB ss Bq. For 
small Kubo numbers, however, the ratio k±/k\\ can be 
very small. A constant value of k±/k\\ was discussed 


before in the context of solar wind observations (see 
Ref. [27] for more details). However, from observations, 
it is not entirely clear what the value of this ratio re¬ 
ally is. It is often assumed that it is very small (e.g., 
k_l/k|| « 0.01 or even smaller) but there are also reports 
that k±/k|| « 1 (see Refs. [28] and [29]). In the current 
paper an explanation for this behavior is provided. 
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